Solitary wave complexes in two-component mixture condensates 
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Axisymmetric three-dimensional solitary waves in uniform two-component mixture Bose-Einstein 
condensates are obtained as solutions of the coupled Gross-Pitaevskii equations with equal intra- 
component but varying intercomponent interaction strengths. Several families of solitary wave 
complexes are found: (1) vortex rings of various radii in each of the components, (2) a vortex ring 
in one component coupled to a rarefaction solitary wave of the other component, (3) two coupled 
rarefaction waves, (4) either a vortex ring or a rarefaction pulse coupled to a localised disturbance of 
a very low momentum. The continuous families of such waves are shown in the momentum-energy 
plane for various values of the interaction strengths and the relative differences between the chemical 
potentials of two components. Solitary wave formation, their stability and solitary wave complexes 
in two- dimensions are discussed. 

PACS numbers: 03.75.Lm, 05.45.-a, 67.40.Vs, 67.57.De 



Solitons and solitary waves represent the essence of 
many nonlinear dynamical processes from motions in 
fluids to energy transfer along biomolecules, as they 
define possible states that can be excited in the sys- 
tem. These are localised disturbances of the uniform 
field that are form-preserving and move with a con- 
stant velocity. The nonlinear Schrodinger (NLS) equa- 
tion iipt + V^V' + tIV'P'/' = is canonical and universal 
equation which is of major importance in continuum me- 
chanics, plasma physics, nonlinear optics and condensed 
matter (where it describes the behaviour of a weakly in- 
teracting Bose gas and known as the Gross-Pitaevskii 
(GP) equation). The reason for its importance and ubiq- 
uity is that it describes the evolution of the envelope 
ip of an almost monochromatic wave in a conservative 
system of weakly nonlinear dispersive waves. Similarly, 
systems of the coupled NLS equations have been used 
to describe motions and interactions of more than one 
wave envelopes in cases when more than one order pa- 
rameter is needed to specify the system. The coupled 
NLS equations have been receiving a lot of attention with 
recent experimental advances in multi-component Bose- 
Einstein condensates (BECs). BECs can excite various 
exotic topological defects and provide a perfect testing 
ground to investigate their physics, because almost all 
parameters of the system can be controlled experimen- 
tally. Topological defects in two-component BECs have 
been predicted theoretically, but there is still no under- 
standing of what the complete families of these defects 
and solitary waves are, nor of their properties, formation 
mechanisms and dynamics. It has also been suggested jlj 
that multi-component BECs offer the simplest tractable 
microscopic models in the proper universality class of cos- 
mological systems and solitary waves in multi-component 
BECs may have their analogs among cosmic strings. The 
goal of this Letter is to find and characterise the fami- 
lies of solitary waves that exist in systems of the coupled 



FIG. 1: (colour online) The dispersion curves of four families 
of the axisymmetric solitary wave solutions of (|3J (A) (red 
line) a = 0.7 and A^ =0.1 (c_ = 0.2738); (B) (black line) 
a = 0.1 and A^ = 0.1 (c_ = 0.6169); (C) (green line) a = 0.5 
and A^ = 0.25 (c_ = 0.3317); (D) (blue line) a = 0.5 and 
A^ =0.1 (c- = 0.3905). The numbers next to the dots give 
the velocity of the corresponding solitary wave. For (D) (blue 
line) these are 0.3, 0.32, 0.34, 0.36, 0.38. AU these solutions 
are VR-VR complexes except for U = 0.58 on (B) branch 
which is VR-RP. The top inset shows the density isosurface 
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complex for a = 0.5, A = 0.25 that is moving with U = 0.3. 
The radii are 6i = 5.194 and 62 = 4.796. The density contour 
plots of this solution are shown in the bottom inset. 
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NLS equation in three dimensions. The implications of 
these solitary waves are wide-ranging, but they will be 
discussed in the context of two-component BECs. 

The simplest example of a multi-component system 
is a mixture of two different species of bosons, for in- 
stance, ''^K-^^Rb 0- Since alkali atoms have spin, it is 
also possible to make mixtures of the same isotope, but 
in different internal spin states, for instance, for ®^Rb 
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0. The multi-component BECs are far from being a 
trivial extension of a one-component BEC and present 
novel and fundamentally different scenarios for their ex- 
citations and ground state The theory for a mixture 
of two different bosonic atoms can be developed similar to 
that for a one-component condensate whose equilibrium 
and dynamical properties can be accurately described by 
the GP equation [Sj for the wave function -0 of the con- 
densate 
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VV(r,t) + K)|V^(r,i)PV'(r,i), (1) 



where m is the mass of the atom, Vb — 'inh^a/m is the 
effective interaction between two particles, and a is the 
scattering length. The GP model has been remarkably 
successful in predicting the condensate shape in an exter- 
nal potential, the dynamics of the expanding condensate 
cloud and the motion of quantised vortices. The family 
of the solitary waves for ^ was numerically obtained in 
1^. In a momentum-energy {p£) plot, the sequence of 
solitary waves has two branches meeting at a cusp where 
p and £ simultaneously assume their minimum values. 
For each p in excess of the minimum pc, two values of £ 
are possible, and £ — > cx) as p — > oo on each branch. On 
the lower (energy) branch the solutions are asymptotic 
to large circular vortex rings. As p and £ decrease from 
infinity on this branch, the solutions begin to lose their 
similarity to vortex rings. Eventually, for a momentum 
Po slightly in excess of Pc, they lose their vorticity, and 
thereafter the solutions may better be described as "rar- 
efaction waves" . The upper branch solutions consist en- 
tirely of these waves and, as p — > oo, they asymptotically 
approach the rational soliton solution of the Kadomtsev- 
Petviashvili Type I equation. 

For two components, described by the wave functions 
tpi and 1/^2, with A^i and N2 particles respectively, the GP 
equations become 
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where mj is the mass of the atom of the ith conden- 
sate, and the coupling constants Vij are proportional 
to scattering lengths via Vij = 2TTh^ aij / rriij , where 
rriij = mirrij / {rrii + rrij) is the reduced mass. 

Two-component one-dimensional BECs have recently 
been considered and various structures have been iden- 
tified such as bound dark-dark, dark-bright, dark- 
antidark, dark-grey, etc. complexes. In higher dimen- 
sions, domain walls and skyrmions (vortons) 9] have 
been identified by numerical simulations. Numerical sim- 
ulations of two-dimensional rotating two-component con- 
densates were performed ^3 ^^"^ structure of vor- 
tex states were investigated. A phase diagram in the 



intercomponent-coupling versus rotation-frequency plane 
revealed rich equilibrium structures such as triangular, 
square and double-core lattices and vortex sheets. These 
simulations give a taste of a rich variety of static and dy- 
namic phenomena in multi-component condensates. One 
would expect the existence of various other families of 
solutions many of which have not yet been detected. 

In what follows I determine the families of three- 
dimensional axisymmetric solitary wave solutions that 
move with a constant velocity U in uniform two- 
component mixture BECs. The trap geometry, relevant 
to experiments, introduces a harmonic-oscillator poten- 
tial in (O together with additional parameters, places 
restrictions on studies of solitary waves and their stabil- 
ity and is irrelevant in the view of our interest to effects 
that occur in large systems. Also, I believe that addi- 
tional physical mechanisms should be introduced only 
after simpler models are well understood. Nevertheless, 
the results I obtain will be relevant to experiments with 
a sufficiently shallow trap, so the linear dimensions of the 
trap are much larger than the healing length. To reduce 
the number of parameters in the system, I will assume 
that the intracomponent scattering lengths and masses of 
individual components in the mixture are equal, so that 
rrii = m and an = a, but the intercomponent scattering 
lengths differ from a. 

To find axisymmetric solitary wave solutions moving 
with velocity U in positive z— direction, I solve 
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where a dimensionless form of lO is used such that the 
distances are measured in units of the correlation (heal- 
ing) length ^ = h/ ^2m^i, the frequencies are mea- 
sured in units 2fii/h and the absolute values of the 
fields IV'ip and |02p are measured in units of particle 
density n = ^i/Vu. Also present in (O are the pa- 
rameter of the intercoupling strength a = V12/V11 and 
the measure of asymmetry between chemical potentials 
= (/zi — /i2)//ii (where we assume that ni > ^2)- 
The dispersion relation between the frequency uj and 
the wave number k of the linear perturbations ((x exp[ik- 
X — iojt]) around homogeneous states, V'ioo, is obtained 
as 



Pk\ (4) 



- k\e + 2^L)) (4c.'' - fc'(fc' +0-200)) 

where P — a^V'ioo'/'foo- The condition of dynamic stabil- 
ity cj(fc) > gives < I If a < — 1, the gas is unstable 
to formation of a denser state containing both compo- 
nents, while if a > 1, then the two components will sep- 
arate. In what follows I will consider < a < 1. The 
values of the wave-functions of the solitary waves at in- 



finity in (PJ are given by 
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FIG. 2; (colour online) The dispersion curves of three fam- 
ilies of the axisymmetric solitary wave solutions of Q with 
a = 0.1 and =0.1. The numbers next to the dots give the 
velocity of the solitary wave solution. The top (black) branch 
corresponds to VR-VR (VR-RP for U = 0.58) complexes. 
The middle (green) branch shows p ys £ for VR-SW com- 
plexes and the bottom (red) branch is the dispersion curve 
of SW-VR (SW-RP for U = 0.58) complexes. The radii of 
the vortex ring solutions are shown on the upper inset as a 
function of U : the two top (red and green) lines give bi and 
62 correspondingly for VR-VR complexes. The two bottom 
(black and blue) lines represent 61 in VR-SW complex and ti2 
in SW-VR complex correspondingly. The bottom inset shows 
3D plots of z)P of the SW-RP complex moving with the 

velocity U = 0.58. 
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and iploa = 1 — mp2oc- the long- wave limit (fc — > 0), 
Q gives two acoustic branches lu± « c±k with the 
corresponding sound velocities c± = ^(V'loo + "^ioo ^ 

Vi^ioo - i^looY + 4aV?ooV'ioo)^^^- The solitary waves 
I seek below are all subsonic, i.e. U < C-. This gives 
a restriction on the asymmetry parameter A^: c_ is real 
only if < 1 - a. 



Each solitary wave complex that belongs to a family 
of the solitary wave solutions for a chosen set of (a, A^) 
will be characterised by its velocity, J7, vortex radii 6^, 
momenta Pi = (0,0, pi), and energy £. The momentum 
(or impulse) of the i— th component is Pi = ^ /[(V'j* ^ 
V'ioo) VV'i — (V'i — '!/'ioo)Vi/'j*] dV. The reasons for replacing 
a more customary defined momentum Pi = ^ / i'l'^'^i — 
tpiVtp* dV with the convergent integrals pi were spelled 
out in for a one component GP equation. Also, similar 
to 0, we form the energy, £, by subtracting the energy 
of an undisturbed system of the same mass for which 
Tpi = const everywhere, from the energy of the system 
with a solitary wave, so that the energy of the system 



becomes 

^ = \j \^i>i? + \^i^2?dv 

+ \j (V-L - Hi?f + (^200 - m^fdV (5) 

+ f j{^^loo-Hl?moo-\i^2?)dV. 

By performing the variation ipi — > ipi -\- Sipi in the inte- 
grals for £ and pi , discarding surface integrals that vanish 
provided 5ipi — s- for |x| — s- cx) and making use of we 
obtain 68 = US{pi +P2), or U — d£/d{pi +P2), where 
the derivative is taken along the solitary wave sequence. 
The same expression is obeyed by the sequences of clas- 
sical vortex rings in an incompressible fluid and by the 
solitary waves of 0. By applying some algebraic ma- 
nipulations that involve some integration by parts, the 
following integral properties can be established: 
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f = 2 y (1 - l^iP - a|V2n(2V'L-V'ioo(V'i + ri))dV 

+ - 1^2 n(2^2'oo- ^200(^^2 + r2))dV{7) 

U{pi+P2) = \ I {I - \iJi\^ - a\i>2\') 

X (3i^L - V-ioo (V-i + V'D - I V-i n dV 

+ 1|(1_«|^^|2_|^^|2) 

X(3^L - V'2oo(V'2 +V'2*) - \^2f)dV. (8) 

These were used as checks on numerical accuracy of the 
solutions I obtain next. 

The axisymmetric solitary waves are found by re- 
writing ||2Jl in cylindrical coordinates (s, 9, z) for the de- 
viations from the solutions at infinity '^i — — V'ioo- 
Stretched variables z' — z and s' — sVl — 2[/2 were 
introduced and the infinite domain was mapped onto 
the box (0, -1) X (— f , f ) using the transformation z = 
tan-i(Dz') and s = ta.n~\Ds'), where D - 0.4 - 0.5. 
Transformed equations © were expressed in second- 
order finite difference form using 100^ grid points, and 
the resulting nonlinear equations were solved by Newton- 
Raphson iteration procedure using banded matrix linear 
solver based on bi-conjugate gradient stabilised iterative 
method with preconditioning. As a — > 0, two compo- 
nents become uncoupled, so the solitary wave sequence 
for each component is following the dispersion curve of 
the one-component GP equation |^ with vortex rings 
(VRs) becoming rarefaction pulses (RPs) as U increases, 
with energy and momentum being appropriately scaled 
by the densities at infinity. For a ^ 0, different com- 
ponents become RP at different critical values of U, so 
variety of complexes become possible. Table 1 gives an 
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example of various transitions from one complex to an- 
other in the system with a = 0.05 and = 0.1. Notice 
that the radii of the vortex rings in VR-VR complexes 
differ with &2 — ^ as f7 — > 0. Also notice, that there 
is a cusp in the dispersion curve E vs pi + P2, since the 
solitary wave complex moving with U — 0.63 belongs to 
the upper branch. 

Table 1. The velocity, U, energy, £, momenta, pi and radii, 6;, 
of the sohtary wave solutions of with a = 0.05 and = 0.1. 
The sequence terminates at (7 = c_ 0.646. 

U £ pi p2 hi 62 complex 

0.55 118 92.3 80.8 1.83 1.50 VR-VR 

0.58 107 81.0 72.6 1.47 0.57 VR-VR 

0.60 101 74.9 69.6 1.12 - VR-RP 

0.63 102 66.2 79.0 - - RP-RP 

As the interaction strength a increases for fixed A^ and 
[/, the radii in VR-VR complexes don't change much, 
although energy and momentum decrease significantly. 
We can compare two such solutions for A^ — 0.1 and 
U = 0.3: the VR-VR complex for a = 0.1 has £ = 387, 
pi = 484, p2 = 410, 61 = 5.197 and 63 = 5.096, whereas 
for a = 0.5 these values are £ = 284, pi = 385, p2 — 271, 
hi — 5.196 and 62 = 5.093. On the other hand, if a 
and U are kept constant and the asymmetry parameter 
A^ increases, the radius of the vortex ring in the second 
component decreases (if A^ = 0.5, a = 0.1, U = 0.3, 
then £ = 284, pi = 503, p2 = 160, 61 = 5.182 and 62 = 
4.363). Fig. n shows the dispersion curves of families of 
the axisymmetric solitary wave complexes for a variety 
of (a,A2). 

In addition to the VR-VR, VR-RP, and RP-RP com- 
plexes that I have just described, two other families of 
the solitary waves came as a surprise. In contrast with 
the VR-VR, VR-RP, and RP-RP complexes, in which 
the solitary wave in each component possesses a large 
momentum, the total momentum of the system in the 
new families is almost entirely belongs to one compo- 
nent with the disturbance of the other component mov- 
ing almost entirely through the coupling with the mobile 
component. I will call this disturbance a "slave wave" 
(SW). This gives rise to two other complexes such as VR- 
SW and RP-SW. Notice that the role of the components 
can be reversed giving SW-VR and SW-RP complexes as 
well. To indicate how the characteristics of these com- 
plexes change with increasing intercomponent interaction 
strength a, the parameters of these solutions are given 
in Table 2 for [/ = 0.3 and A^ = 0.1. FigH shows the 
dispersion curves of all three families of the axisymmetric 
solitary wave solutions for a = 0.1 and A^ = 0.1. 

Table 2. The intercomponent interaction strength, a, energy, 
E, momenta, pi and radii, hi, of the solitary wave solutions of l3l 
with U = 0.3 and = 0.1. The sequence terminates (c_ = U) at 
a ^ 0.65. 



a £ pi p2 hi 62 complex 
0.2 178 414 1.52 4.977 - VR-SW 
0.4 142 330 7.76 4.673 - VR-SW 
0.6 119 268 35.1 4.122 - VR-SW 

a £ Pi P2 hi 62 complex 

0.2 140 1.12 330 - 4.787 SW-VR 

0.4 96.7 4.89 237 - 4.354 SW-VR 

0.6 63.0 17.5 174 - 2.986 SW-VR 

Finally, when the solitary waves are found, there is 
a need to elucidate the following topics. Stability. The 
Derrick-type argument used in can be applied together 
with the integral identities I©-®, to indicate the sta- 
bility of solitary waves as long as dU/d{pi + P2) < 
along the family. This condition is satisfied for all solitary 
waves discussed here except for the last entry of Table 1 . 
Solitary waves in 2D. Similarly to 3D complexes, solitary 
waves in 2D consist of the combinations of a pair of two 
vortices with opposite circulation (VP), 2D rarefaction 
pulse, or a "slave wave," giving rise to VP(RP)-VP(RP) 
and SW-VP(RP) complexes. Formation. The mecha- 
nisms of solitary wave nucleation in two-component con- 
densates are similar to those in one-component conden- 
sates 1^3- III particular, the moving objects (ions, laser 
beams etc) shed VR(RP)-VR(RP) complexes in 3D and 
VP(RP)-VP(RP) complexes in 2D, when the local speed 
of sound is reached somewhere on the surface of the ob- 
ject. I will ellaborate on these issues in details elsewhere. 
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